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Abstract 

The N-quantum approach (NQA) to quantum field theory uses the complete and 
irreducible set of in or out fields, including in or out fields for bound states, as stan- 
dard building blocks to construct solutions to quantum field theories. In particular, 
introducing in (or out) fields for the bound states allows a new way to calculate 
energy levels and wave functions for the bound states that is both covariant and ef- 
fectively 3-dimensional. This method is independent of the Bethe-Salpeter equation. 
In contrast to the Bethe-Salpeter equation, all solutions of the NQA are normalizable 
and correspond to physical bound states. In this paper we use the NQA in one-loop 
approximation to calculate states of the relativistic hydrogen atom and analogous two- 
body systems to illustrate how our new method works. With additional terms in the 
in field expansion we find systematic corrections beyond the Coulomb interaction. 

1 Introduction 

The tour de force experiment of Pohl, et al. jl] and [2] provided new motivation for us to 
continue work on the NQA of calculating bound state properties. The discrepancy between 
the proton structure measured in hydrogen and that measured in muonic hydrogen has 
three possible causes, (1) new physics, (2) inadequate QED calculations, or (3) incorrect 
description of the interaction with the proton. To study (2), we are developing a new way 
to do the QED calculations. In this paper we introduce our new method of calculation, 
but do not carry the work to the level necessary to resolve the discrepancy found in the 
Pohl, et al, and Antognini, et al, experiments. Other attempts to resolve the discrepancy 
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can found from the citations to these articles. We cite one such calculation [3j that also 
examines the second case. In this paper we present a one-loop calculation of the energy 
levels and wave functions of ordinary hydrogen and muonic hydrogen. The NQA, based on 
Haag's expansion [4j of interacting fields in terms of asymptotic fields, also can be applied 
to other two-body systems, such as the (ejl) and (fifl) systems. 

Haag's original expansion did not take account of bound states. We added in (or out) 
fields for each of the bound states that we take as stable in our approximate treatment. For 
the case of the hydrogen atom we add an in (or out) field for every state of the hydrogen 
atom. 

2 Goals of this paper 

Our main goal in this paper is to develop a new method of calculation of the energy 
levels and wave functions of relativistic bound states. This method has been described 
previously j5j [61 [8] , but it has not been developed sufficiently to account for high-order 
radiative and recoil effects that are relevant to the analysis of high-precision spectroscopic 
measurements, such as those in ordinary and muonic hydrogen. Among the advantages of 
this method for (2-body) bound states such as the hydrogen atom are the introduction of 
both masses as independent parameters, rather than via the reduced mass. Since we are 
interested in muonic hydrogen as well as ordinary hydrogen, taking the proton mass as 
an independent parameter is important. We also solve integral equations that incorporate 
radiative and recoil effects without perturbation theory. References to other applications 
of the NQA are in Of particular relevance to our present paper is |6j, which is the first 
paper to give the spectator equation for a two-body bound state. There is related work by 
Kallen, |10] . Yang and Feldman [llj and, for bound states, by Gross [12j. 

In this paper we present a one-loop approximation for the relativistic bound state 
that reduces to the Dirac equation for large proton mass. As stated above, we introduce the 
proton mass as an independent parameter. This is relevant in the case of muonic hydrogen 
where the mass ratio m^/M p ~ 1/10 rather than the ratio m e /M p ~ 1/2000 for ordinary 
hydrogen. We use the Coulomb potential as the binding mechanism and ignore magnetic 
interactions and renormalization effects in this paper. We do not calculate energy levels in 
high accuracy in the present paper. Rather, we describe our new method of calculation. 
In later papers, in addition to our systematic development of the NQA, we will give a 
unified calculation that includes all quantum electrodynamic corrections to the hydrogen 
spectrum up to a given order, rather than adding various corrections piece by piece. In our 
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next paper we will pay particular attention to any differences between our results and the 
usual calculations to see whether this new method resolves the muonic hydrogen anomaly 
concerning the proton charge radius as inferred from measurements of the Lamb shift in 
ordinary and muonic hydrogen. We will also find a set of coupled integral equations that 
include all terms up to a relevant order, rather that adding corrections term by term. 

3 Asymptotic fields and the Haag expansion 

The in (out) fields have free field commutators, obey free equations of motion, and the 
different in (out) fields commute or anticommute with each other everywhere in spacetime. 
Each of these sets of asymptotic fields by themselves is completely known once the masses, 
spins, and quantum numbers of the fields in a given set are given. Thus either set serves 
as a collection of standard building blocks to construct solutions of the operator equations 
of motion. For the present paper we ignore the difficulty that the asymptotic limits for 
the charged fields do not exist. In a recent paper [13J we found modified charged fields for 
which the asymptotic limits do exist; however we do not use the modified charged fields 
here. 

To solve the equations of motion using the Haag expansion, expand the fields that 
appear in the Hamiltonian or Lagrangian in normal-ordered series of in (out) fields. To 
determine the c-number amplitudes (the Haag amplitudes) that are the coefficients of the 
normal-ordered terms, insert the expansion in the operator equations of motion, renormal- 
order, and equate the coefficients of corresponding (linearly independent) normal-ordered 
terms. The relevant Haag amplitudes are the wave functions of the bound states. 

The method based on the Haag expansion is entirely independent of the Bethe- 
Salpeter equation p3]. In contrast to the Bethe-Salpeter approach, in the NQA there are 
no spurious solutions and no negative norm amplitudes. The NQA can be used for bound 
states in relativistic theories even though the amplitudes depend only on the same number 
of kinematic variables as nonrelativistic wave functions. In particular, there are no relative 
times in the relativistic version of the NQA. With all terms in the in field expansions 
allowed by conservation laws, the Haag expansion should be equivalent to the interacting 
field theory. This results in an infinite set of coupled equations. To get a tractable set 
of equations, we terminate the Haag expansions, keeping a finite set of terms for each 
interacting field. For quantum electrodynamics, the case relevant here, the smallness of 
a provides a rationale to terminate the series; each vertex has a factor of \fa. For this 
paper we restrict to the simplest terms in the Haag expansion that give an equation for 
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the bound state. 

We need not construct the in field for the bound state. We assume possible bound 
states and introduce in fields, characterized by their mass, spin, and other quantum num- 
bers, for each. The equations of motion for the interacting fields give equations for the 
bound state amplitudes; if there is a solution for a given bound state amplitude, then 
the corresponding bound state exists. We take "bound state amplitude" as a synonym for 
"wave function." 



4 Relativistic model of the hydrogen atom 

The fundamental fields are the electron, e a (x), muon, fj, a (x), proton, p a (x), and photon 
vector potential, A^{x) fields. These fields obey the operator equations of motion (for this 
paper we drop renormalization counter terms), 

(t^-m)e(x) = |[4(x),e(x)] +) (1) 

(i 0- m^)u{x) = -yL(x),(i(x)]+, (2) 

(i@-M)p(x) = ~[fi(x), P (x)] + , (3) 

d»d ■ A - 8 ■ 8A» = |([e(x), 7 M e(x)]_ + \p,(x),-fn{x)}- - \p(z),-fp(x)]-) (4) 

where Eq.Q follows from 

d v F^{x) = |([e(x),y*e(x)]_ + [/x(x), 7 ^(x)]_ - \p(x),-fp(x)].) (5) 

and F^ u (x) = d^A u {x) — d^A v {x). We choose the masses of the electron, muon, proton, 
and hydrogen atom in states i as m, run, M, and Mj, respectively. The equations are 
symmetric under e O p, too ma- 

We use the Haag expansion to expand the interacting fields appearing in the equa- 
tions of motion in terms of in fields. We truncate the series, keeping the first term involving 
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the hydrogen bound state in fields, hf 1 , 




x — y, x — z)i 




(6) 



e( x ) = eW(x) + d "y^ z ■ hP\z)i A fph(x -y,x- z)p^ n \y) : (7) 

i 

p{x)=p^\x) + Y J I d 3 yd 3 zf- eh (x - y,x - z) : e^\y)i A h W( z ) : ( 8 ) 

i 

p(x) = pW(x) d 3 yd 3 z : hf n) \z)i A e (in) {y) . /gfe(x _ y , x _ Z ) ( 9 ) 

A"(x) = ^"^(x) + f d 3 yd 3 z[:p^(y)f^(x-y,x-z)p (m \z): 



+ :e^\y)f^ e (x-y,x-z)e^\z) :] 



(10) 



where the J2i i s the sum over the various hydrogen states which, for simplicity, we took as 
scalar, 



and we label each amplitude by the in fields in each term in the expansion of the interacting 
fields, and keep this label for the terms in the adjoints of the interacting fields. 

We used translation invariance to write these forms of the expansions. Lorentz 
covariance gives the transformation properties of the amplitudes: 



We choose spectroscopic notation for the states of the hydrogen atom that is adapted to 
treating the proton spin on the same basis as the electron spin. We use F, L, S for the 
total angular momentum (an exact quantum number), the orbital angular momentum, 
and the lepton-proton spin, respectively. With the principal quantum number, n, we label 
states as nLg. (Our choice differs from the usual choice that couples the orbital angular 
momentum, L, to the electron spin, S e , then couples J = L + S e to the proton spin to get 
F, and labels the states as nLj.) 



f ph (x,y) = 1 °f ph (x,y) 7 
hh(x,y) = J°fl h (x,yh°. 



(11) 
(12) 



S(A)f ph (x,y)S(A)- 1 = f ph (S(A)x,S(A)y). 



(13) 
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Because the contractions, (0|e( m )(a:)e( m )(y)|0), etc. are simpler in momentum space 
than in position space, we continue our analysis in momentum space. To go into momentum 
space, we use 

e(x) = J d 4 qe(q) exp(— iq ■ x) (14) 

and analogous formulas for the other fields. We leave tildes off the Fourier-transformed 
fields. The equations in momentum space are 

(4- m )e(q)= e - j d 4 k[fi(k),e(q-k)} + (15) 
W-™>(</) = | J d i k[4(k),f i (q-k)] + (16) 

^-M)p( P ) = - e - J d%m,p(p-m+ (17) 

-Wk ■ A{k) + k 2 A»(k) = | J d 4 q'([e(q'),^e(k - q')\ (18) 

+ W),-fM - - W)^p{k - </)]) (i9) 

To avoid subscripts we use h for ordinary hydrogen (electronic hydrogen) and H for muonic 
hydrogen. For this one-loop approximation we choose the Coulomb gauge. We expand the 
interacting fields in normal-ordered products of in fields. For the one-loop approximation 
to the amplitude in which e ~: p^ 1 

n) h (in) 

:, we keep terms with up to three in fields in the 
Haag expansions for e and A and one contraction. These terms are : pp^Tph :, : A :: Aph :, 

pp~^TpAh :, where the overbraces stand for contractions. The term from : hh :: hp : is 
much higher order because the hydrogen atom has zero charge. The Haag expansion for 
the electron field is 

e(q) = e m (q) + £ J d 4 pd 4 b5{p + Q ~ b)fp hj (p, b) : p in (p)hf(b) : 
j 

+ j d 4 pd 4 b5(p + q + k- b)f Aph . (p, k, b) : A^-k)p m (p)hf(b) : (20) 
e(q)=e in (q) + J2 j d 4 pd 4 b5(p + q - b) : h m \b)p m {p) : f ph .(p,b) : 

3 

J d 4 pd 4 b5(p + q + k-b): h in \b)p in (p)A^k) : f Aphj (p, b) :, (21) 



+ 



where f ph .(p,b) = J 0T f ph . (p, b)j°. 
We chose this parameterization so that 

(e(?)p(p)|0), hl(b)\0)) = 5(g + p - b)f phi (p, b)(p + M)e( P °)5( P 2 - M 2 )9(b°)5(b 2 



Mi 



(22) 



There are analogous expressions for the muon and proton fields. For the photon field, 
A »(k) = A» in (k)+ [ d 4 pd A p'5(k-p-p')[:p m (p)f^(p,p')p m (p') : 



- ■■ e m {p)r eel {P,p')e m {p') : - : tf n {p)Q{p,p')n m {p') :]■ (23) 

After re-normal-ordering, we find the one-loop equations for the two main amplitudes for 
any state of the hydrogen atom, 



m)f e (p, b) 



2(2vr) 3 J 2E p , ' (p - p 



(l/-</-M)f p ( q ,b) 



2(2vr) 3 J 2e p/ 1 (b-p' - p) 
1 f W ufptf,b), , T 



2(2tt) j ^ 

e 2 [d s q> u Uq',bf 



2(2tt) 3 J 2E q > 1 (b-q'-q) 2 



(24) 



(25) 



where f e (p,b) = f ph (p, 6)(j/ + M) T , f p (q,b) = fehMW + ™f , E p = ^P 2 + M 2 , 
e q = \/<i 2 + m? , p is the energy-momentum of the on-shell proton, and q = b — p is 
the energy-momentum of the off-shell electron. Note that, by construction, f e (p,b) obeys 
the subsidiary condition f e (p,b)(p / — M) T = and f p (q,b) obeys f p (q,b)(q / — m) T = 0. 
Unlike the Bethe-Salpeter approach, we have arrived at a pair of coupled equations that 
describe the bound state. They are explicitly symmetric under subscript e -H- p and mass 
m O M interchange. These two equations differ from those found in |15| and |16| where 
Bethe-Salpeter equations with one on-shell particle are found, but we will show that in 
a certain approximation they reduce to their Bethe-Salpeter counterparts. As far as we 



know, the exact properties of Eqs. (24) and (25) are unexplored. 



The terms on the right hand side of Eqs. ( 24 ) and ( 25 ) are expressed in diagrammatic 



form in figures la and lb Heavy lines are off-shell and light lines are on-shell. Point vertices 
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(a) 



(b) 



Figure 1: Graphs for the right hand side of the electron equation of motion. Heavy lines are 
off-shell and light lines are on-shell. The dashed line represents the bound state (hydrogen atom). 
The empty circle represents the amplitude f e in (a) and f p in (b). The left fermion line is the 
electron and the right line is the proton. Similar graphs exist for the proton equation. 

represent the substitution of an off-shell interacting field in terms of other interacting fields 
via the relevant equation of motion. These are the fundamental QED vertices. Circles 
indicate the use of the Haag expansion to express off-shell interacting fields in terms of in 
fields with a Haag amplitude coefficient. We will use these diagrams in section 9 to show 
how higher order corrections are calculated. 

In a covariant gauge these equations are explicitly Lorentz covariant for an arbitrary 
state of motion of the hydrogen atom. The general expansions have the support of the in 
fields on both mass shells. Here we keep only the mass shell that does not lead to extra 
pairs of particles. 

5 Normalization of the wave functions 

The asymptotic fields diagonalize conserved observables such as the Hamiltonian, the mo- 
mentum operators and various charges. We can represent any conserved quantity O in 
terms of either the interacting fields or the asymptotic fields (either the in or out fields), 



where £j in an interacting field and "ipf 1 is an asymptotic field, including the in fields 
for bound states. We do not include weak interactions in our analysis of the hydrogen 
atom; thus both the number of electrons and the number of muons are conserved. We use 
the number of electrons, N e , to find the normalization condition for the hydrogen wave 




(26) 
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function, fph(p, b), in which the electron is off-shell. The only interacting field that carries 
electron number is the electron field, 

N e = f e^(x)e(x)d 3 x. (27) 



The contribution to the electron number from the hydrogen atom in a given state comes 
from the terms in N e that are bilinear in the hydrogen atom in field in that state. From 
the in field expansion of e(x) we find the orthonormalization condition 

J M ^Th:[f pH . l (p, b')^f pH] (p, b)] = 5 fj 2E b 5(h' - b). (28) 

6 Interchange of the on-shell and off-shell particles 

The equal-time anticommutators relate the Haag amplitudes with the lepton off-shell to 
those with the proton off-shell. These relations follow from the vanishing of the coeffi- 
cients of each (linearly independent) normal-ordered product of in fields in the equal-time 
anticommutators. Most of the relations involve Haag amplitudes for terms with higher 
degree normal-ordered products than we have considered here; however for the equal-time 
anticommutator [e,p]+ = there is an approximate relation that involves only terms that 
we considered here, 

[feMX^fW + [fpH,(vM^^) T U = o, (29) 

with the constraint p + q = b. Thus, the Haag amplitude with the lepton off-shell is 
simply related to that with the proton off-shell. The two amplitudes determine each other 
uniquely. 



Using Eq. (29), we can simplify Eq. (24) to 



(P-/- -»)/.<»») = (4 / Mf , (30) 

which matches the Bethe-Salpeter equations of [15] and |16J . For the present paper, we 
choose the hydrogen atom at rest, b = (M^, 0), which explicitly breaks Lorentz covariance 
to rotation covariance. Keeping the main mass shell, and dropping the magnetic interaction 
terms, we have 

(i°Mh -i/-m)/ e (p) = -/^ / ffl°r^%{l°) T ^+M) T , (31) 

(27r) d J Ihpi |p-p'r 
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where f e (p) = f e {p; M Hi ,0). 



7 Solution to bound state equation 



The purpose of the following sections is to find a method for solving Eq. (31 ) which can be 



extended to solve Eqs. (24) and (25). For the sake of simplicity in this work, we will put off 
solving our more complicated coupled equations in a future paper. We acknowledge that we 
are solving an equation that has already been studied extensively in the literature, but our 
purpose is to develop the NQA framework for high precision calculations. We therefore 
take an approach that differs from the typical perturbative method. The methods for 
finding higher order corrections will be discussed in section 9. 

In this section, we focus on binding due to the Coulomb interaction. We chose 
Coulomb gauge to simplify our calculations and to allow comparison with the usual solution 
of the Dirac equation for the hydrogen atom. To keep the notation general for any two- 
particle system, we label the constituents m\ and mi and the bound state m& in this 
section. We solve the bound-state equation 

(7 mfo _y_ mi)/e(p) = j^f ^y7 ^(p,p , )/ e (p , )(7°) T (y+^2) T (32) 

2-Ey 

where m^ = m\ + mi + e«, m^ is the energy of the hydrogen state, e« < is the binding 
energy of the atom, V(p,p') = — l/|p — p'| 2 , and Ejfi = \Jp 2 + mf . A similar equation 
is solved in [TJ [8] in the non-relativistic limit. We solve the equation numerically without 
taking a non-relativistic limit. 

We can think of mi as the mass of the proton and m\ as the mass of the lepton, but 



the N-quantum equations that give Eq. (31) are symmetric under m\ o m-i together with 
e o — e and our calculations reflect this. 

Before solving this equation, we show that it reduces to the expected Dirac equation 
in the large mi limit. The factor mi) T /2E p > — > (1 + 7°)/2 in the potential in Eq. (31) 



reduces the 4x4 system of equations to a 4 x 2 system with the usual Coulomb potential. 



From Eq.(31 ), using q = b — p, we find 

(£/ _ 7 . q _ mi _ 7 <V )/e = Q) (33) 

where V is the Coulomb potential. Because this equation comes from a covariant formu- 
lation, we have to multiply from the left by 7 = j3 to get the usual form of the Dirac 
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equation for the hydrogen atom, 



( a -q + m 1 (3 + V)f e = Ef e , (34) 

using 7° = /3 and 7°7 1 = a 1 . 
7.1 Form of the matrix 

To solve this equation, we break the 4x4 matrix down into four 2x2 matrices: 

/.(p) _ ( • 4 < P > B(P » ) . (35) 
\C( P ) D(p)) 

Next we introduce the partial wave expansion of the operators and the amplitude. Each of 
these 2x2 matrices can be written as a product of a spin- angle part and a radial function. 
For example, for a specific eigenstate we can write 

A(p) = & L F s mF (n)g L ( P ), (36) 

where W[™ F (Q) is the spin- angle function, #l(p) is a radial function and p = |p|. The 
most general solution is a sum over all possible eigenstates. The spin-angle function is 
given by 

^ s mF (n) = ^< LS;m L m F -m L \Fm F > Smji ._ mi >WM)- (37) 

where 4>g ms is the total spin state of the constituents, YL mL (9, 4>) is a spherical harmonic, 
and < LS;m,LmF — m^FrriF > is a Clebsch-Gordan coefficient. The spin state can be 
either a singlet or a triplet. We express these in terms of 2-component Pauli spinors: 



0oo = ^ Wt) ® XU) - m ® X(t)) = ^= ( _° x M (38) 



0ii = V'(t) ® x(t) 



1 




(39) 







l-i 



Hi) ® x(4) 




1 



(40) 
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0io = ^ Wt) ® X(i) + ® x(t)) = ^ ( ° J j ■ ( 41 ) 

We expect our matrix wave function to be analogous to the direct product of an 
electron and a proton spinor, 



* s *e®*j=i ^ e0 ^ T ^®("-p^) T y (42 ) 

(<T • pV>e) <8> V'p ( CT • P^e) 8) (<T • p^p) T ' 



where \I/ e and ^ p are free Dirac spinors for the electron and proton respectively, and ip e 
and tp p are their upper components. Our wave function must also satisfy the auxiliary 
condition 

/ e (p)(j/-m 2 ) T = 0. (43) 
With these two things in mind, we use the form 

f(v)=( ^s mF 9L{p)l S(pW L F s m -^-p) T g L (p) \ (M) 

Je{V) { tr-wffi'hM s(p)*-p^(*-p) T h L (p) )' 

(2) 

where S(p) = pj {E p -\-rn2), and p is the unit vector in the direction of p. We constructed 



this wave function to satisfy Eq. (43). This wave function is also a parity eigenstate, 



7°/ e (-p)7 0T = (-l) L /e(p). (45) 
7.2 The coupled radial integral equations 



Using Eq.(44) and the LHS of Eq.(31) gives 



where 



LHS= L r U L r a (46) 

\ -^21 -^22 / 



L n = [(m b - - m 1 )g L {p) + ph L {p)]&£™ F 
L12 = S(p)[(m b - E® - mi)g L (p) + ph L {p)]&[™ F {* • p) T 
L21 = -\pg L (p) + (m b + mi - E^)h L {p)]a ■ f>&[™* 
L22 = S(p)[pg L (p) + (m b + mi- E®)h L (p)]v ■ p^ F {<r ■ pf 
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The RHS becomes 



RHS = ( Ru Rl2 ) (47) 



where 



Ru = J dV^|^[(4 2) + ^W'lT + S( P ')p^ F L r(a ■ p') T (cr • p) T ]g L (p') 

R12 = [ d 3 p'^^S{p)[(4 2) + ™ 2 )<* • P) T + S(p')p&' F L r{° ■ P'fMp') 
J 2E ( p V 

Rn = ~f d 3 p ,Vij ^f[(E^ + m 2 )a ■ 

+ S{p')pa ■ v'^'ls F ^ ■ P')V • P) T ]h L ( P ') 
R22 = ~ I d 3 p'^^S(p)[(E^ + m 2 )a ■ p'<r («t • P) T 



where ^ ^5 F = <3f£™ F (Q,'). At this point, there is an apparent redundancy in the four 
equations. Right multiplying the upper right and lower right component equations by <r-p 
and dividing by S(p) results in the upper left and lower left component equations. Since 
we reduced the number of independent radial functions in our matrix to two by demanding 



+ S(p')pcr-p'^' F L ^(a.p') T ]h L (p' 



that it satisfy the auxiliary condition, Eq.(43), we expected this redundancy. We will focus 



only on the left components for the remainder of this discussion. 

To keep this analysis general, we must find the action of the a ■ p operators on the 
spin- angle functions, 



- • = £ CF LSL%^VS' F (48) 

L'S' 

• pf = £ C^lTvs^T (49) 



where C F ™ L F , S , are coefficients that can be determined explicitly and tabulated, a ■ p is a 
pseudo-scalar operator and must change L by ±1, i.e. \L — L'\ = 1. Other properties of 
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these coefficients are 



n Fm F _ n T Fm F 
^ LOUS' — ~° LOL'S'i 

n Fm F _ n T Fm F 
°L1L'0 — ~° L1L'0> 

n Fm F _ n T Fm F 
°L1L'1 - ° LlL'li 

^Fmp ^Fmp 

^LSL'S' ~ ^L'S'LS' 



G L5L'5' G L'S'L"5" = d LL"OSS", 

US' 

E nT Fm F „ T Fm F 
U LSL'S ,L - L'S'U'S" - d LL"0sS"- 



US' 



These properties are useful when using our general equations to determine specific 
The partial- wave expansion of the potential is, 

^ oo 

V(p,p') = ^ + 1)Vl(p,p')Pl(cos8 pp ,) 

71 L=0 

oo L 

= -EE v L (p, P ')Yi mL (n')Y LmL (n). 

L=0 m,£,=— L 
Using the orthogonality conditions, 



J dxP L ,(x)P L (x) = l^pfiLV 



J mi J ra T i 



J dnYZ mL (n)Y L , m , L (n) = 6 LL ,s n 

the components of the partial wave expansion in terms of the potential are, 

V L (P,P) =tt 2 J^dxP L (x)V(p,p'), 

where x = cos Q vv > . The orthogonality relation of the spin- angle functions, 

J dQTT[^ L F s m ^^]=hL',S S s' 

are also useful. 



14 



Using Eq. (51) and Eq.(52) we find the left components on the right hand side 



Rn = IJ ^[(EP+m 2 )V L (p, P 'W L F s n 

+ S{P')P £ £ ^L^T^VvM^S^W) (55) 



L 

L'S' L"S" 



R21 = ~l \ ^[(4 2) +- 2 )E^^(p,p')« j 

J 2E W us , 



nip 

jl\^p ■ T- '"-2; ^ ^LSL'S'VL'yfiV }^L>, 
Frnp n T Fm F n T Fm t 

L'S' L"S" L"'S"' 



(56) 



We remove the spin-angle functions by multiplying the top left by (^™ J )^ and the bottom 
left by (<r • p£^™ J )t, taking a trace, and integrating over using Eq,(|54[), The resulting 
equations are 



(m b - - mi)g L (p) + ph L (p) = 

2 / ^H 2) +^) v ^p') + s(p')p^( cT l^ ( 57 ) 

^ J 2E P ' L'S' 

- P9l{p) - (m b + m 1 - Ef>)h L ip) = 



1 1 7^[(4 2) +- 2 )E(^) 2 ^ ( p,p') 

J 2,E pl us , 



+ S(p')p^2 E E ^LSL/S'^ 7 L'S^V'S"^ 7 L''S^'L'''S'''^LSL'''S'''^L''{P,p')]hL(p')- 
L'S' L"S" L"'S"' 

(58) 
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Substituting m& = e + mi + m<i gives 
= (^p 2) - rn 2 )g L (p) - ph L (p) 

+ 2 / ^[(4 2) +^)vb(p>i0+^)p£^ 

(59) 

e^i(p) = (^p 2) - m 2 - 2m 1 )h L (p) - pg L (p) 

+ l I ^[(4 2) + - 2 )E(^) 2 ^(p,pO 



i'p' 2 

L'S' 



L'S' L"S" L"'S"' 

(60) 



7.3 Specific cases of the bound state equation 

As shown earlier, our equation reduces to the Dirac Coulomb equation in the large-7712 
limit. Here we show this reduction for each partial wave. The last term in both equations 
goes to zero and the equation simplifies to 

£9l{p) = -phL{p) + J dp'p' 2 v L (p,p')g L (p') (61) 

eh L (p) = -2mih L (p) - pg L (p) + / dp'p' 2 ^2(C^ s ,) 2 v L >(p,p)h L (p) (62) 

J US' 

where vl{p,p') = ^Vl(p,p')- Again, we find the momentum space Dirac equation for an 
electron moving in a Coulomb potential. 



We can use Eq.(50) along with some general properties of the coefficients to simplify 



our equations in some specific cases. For the case where S = 0, S' must be 1, and we 



can use Eq.(50) to greatly simplify the sums in the last term of the second equation. The 



16 



result is 

£9l(p) = {Ef ] - m 2 )g L (p) - ph L (p) 

+ - f ^[(4 2) +^)^(^pO + 5(pOpE( CT ^i) 2 ^(p,p')]5l(p / ) (63) 
eh L {p) = (E^ -m 2 - 2m 1 )h L (p) - pgi(p) 

+ H ^CH 2) + «2) Y.( C l™Ji) 2 Vl>(p,p) + S(p')pV L (p, P ')]h L (p'). (64) 
J 2E p> v 

For S = 1, L = J, we know L' = J ± 1 and S" = 1. Our simplified equations are 
e 5L(p) = {Ef ] - m 2 )g L (p) - ph L (p) 

„/2 



7 2E p> V 
eh L (p) = (E^ -m 2 - 2m 1 )h L (p) - pg L (p) 

+ 11 ^[(E^ + m 2 )J2(ClrJi) 2 VLM (66) 

* J 2Ey vs , 

Finally, we have the case where 5=1 and L = J — 1. In this case L' must be equal to J, 
and the remaining sum of the squared coefficients over S' is 1. The simplified equations 
are 

£9l(p) = {E^ - m 2 )g L (p) - ph L (p) 

+ ~ I ^C((4 2) + m 2 )V L (p,p') + S(p')pVj(p,p')}g L (p') (67) 
v 

eh L (p) = (E^ -m 2 - 2m 1 )h L (p) - pg L {p) 

P' 

+ S(P')P E E E Cl™f s ,C TF jZ^^ (68) 



S' L" S'" 
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Note that even without the inclusion of a hyperfine spin-spin coupling term there is 
a difference between the nS® and the nS\ equations. The former's state equations, found 



from Eq.(63) and Eq. (64), are 



tgo(p) = (£p 2) - rn 2 )g (p) - ph {p) 



+ - / ^y[(4 2) +™ 2 )V (p,p') + 5(p , )p^i(p,P , )]5o(p') 



(69) 



eMf) = (^p 2) - m 2 - 2m 1 )h (p) - pg {p) 



+ ~ l +m 2 )V 1 (p,p') + S(p')pV (p,p')]h (p') 



(70) 



and the latter's, found from Eq. (67) and Eq.(68), are 
tgo(p) = (M 2) - m 2 )g (p) - ph (p) 



+ 



2 f dp'p' 



(E^ + m 2 )V (p,p') + S(p)pVi(p,p)}g (p') 



(71) 



^o(p) = {E^ -m 2 - 2mi)h (p) - pg {p) 



+ 



2 f dp'p' 



{E^ + m 2 )V 1 (p, P ') + S{p')p{ l -V {p,p') + lv 2 (p, P '))]h Q (p'). 



(72) 



Because p ~ afx, where \x is the reduced mass, the terms containing S(p') are very small. 
For this reason the splitting between the energy levels of these two states created by the 
dissimilarity in the equations is very small. For large m 2 , the potential terms with S(p') 
are smaller by a factor that is 0(a 2 (mi/m 2 ) 2 ). 

7.4 Corrections to the approximation using the reduced mass 

The NQA introduces both the light and heavy particle masses independently, rather than 
introducing the heavy-particle mass via the reduced mass. In the small-p approximation, 
the NQA coincides with the reduced mass approximation, but for larger momenta the 
reduced-mass approximation fails. We show this by examining the kinetic terms in the 
bound-state equations, 
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E 



(2) 



m 2 



-P 



n(2) 



The eigenvalue of this equation is 



-p 




a 


= e 


5 


m.2 — 2mi 




h 




h 



(73) 



e = a/F 2 + m 2 ~ m 2 + VP 2 + — ?™1 



^ — r-P • 

2/x 8m^?n,2 



For the Dirac equation the kinetic terms are 



-p 




9 


= e 


5 


_ -P -2/i _ 




h 




h 



(74) 



The eigenvalue for the Dirac equation is 



y/p 2 + jj? - [i 



P P 
2/i 8/i 2 ' 



where /u = m\mil (jn\ + 7712) is the reduced mass. (This result was found earlier by 
Raychaudhuri. [S]) We note that there are further mass dependencies in the potential 
term, but do not discuss them here. 

8 Numerical results 



We briefly discuss some of our numerical results here. The purpose of this section is to 
show that our procedure and numerical calculations yield results that are consistent with 



standard calculations. We are aware that solutions to Eq. (30) are already well known, 
and we merely intend to show that our procedure does not return any erroneous results. 

We solved the integral equation numerically for several states. We used a grid with 
1200 points per equation and converted the integral eigenvalue equations into matrix eigen- 
value equations. We handled the singularities at p = p' in the kernels with Lande sub- 
tractions |17j . We excluded momenta close to infinity to avoid infinities in our discretized 
integral equation. The wave functions are extremely close to zero well before our cutoff 
is imposed. We made our equations dimensionless by dividing by m\ and expressed the 
coupled equation in terms of the dimensionless parameter £ = milmx. Our results agree 
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n 


NQA 


Dirac 


1 


-13.59847 


-13.59847 


2 


-3.39963 


-3.39963 


3 


-1.51094 


-1.51094 



n 


NQA 


Dirac 


1 


-2.528506 


-2.528527 


2 


-0.632130 


-0.632134 


3 


-0.280946 


-0.280947 



(a) Electronic hydrogen uSq 



(b) Muonic hydrogen uSq 



Table 1: Energy eigenvalues for electronic and muonic hydrogen states. The table on the left gives 
electronic hydrogen levels in units of eV and the right table gives muonic hydrogen levels in units 
of keV. We give Dirac eigenvalues for the Dirac-Coulomb equation with the reduced mass. We 
found NQA values numerically from the the NQA integral equations with a Coulomb potential. 



with those of the Dirac-Coulomb equation with the reduced mass. With a higher precision 
we expect our results to differ from the Dirac-Coulomb equation, because our equation 
contains effects of the proton spin that are not found in the Dirac equation. 

We found a rough estimate of our uncertainty by finding the eigenvalues with 800, 
1000, and 1200 grid points and analyzing the stability of the eigenvalues. We conservatively 
estimated our uncertainty to be 0.01 meV for electronic hydrogen and 2 meV for muonic 
hydrogen. 

We give comparisons of the NQA electronic hydrogen eigenvalues and Dirac-Coulomb 



eigenvalues for the tiSq states in table |l| [(a) These values are nearly identical and the re- 
sults indicate that we may have overestimated our uncertainty. We give the same compar- 



isons for muonic hydrogen in table 1 (b) These values are similar, but differ significantly 
for the lower eigenvalues. The NQA energies in the ground- and next lowest-states are 
higher than the Dirac energies by 21 meV and 4 meV, respectively. It is possible that our 
numerical calculations failed for these two particular eigenvalues, or we may have under- 
estimated the uncertainty. We plan on achieving a higher degree of precision in the future 
to investigate such concerns. 

Using the same method of estimating the uncertainty as before, we conservatively 
take our uncertainty to be 0.01 meV. As in the case of muonic hydrogen, there are some 
discrepancies in the first two eigenvalues. The first and second values are larger than the 
Dirac energies by 0.16 and .03 meV respectively. The higher eigenvalues are consistent 
with the Dirac energies. 

We also calculated the energies of the nSl states. They are identical to those shown 
in tables l| |(a) and (b) for some nS® states. We need higher precision to study the energy 
splitting in these states caused by the differences in the NQA equations. 
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0.010 0.015 



(a) ISq state wave functions 



0.020 0.025 0.030 i 




0.010 0.015 0.020 0.025 0.030 i 



(b) 15*0 wave function differences 




j 1 ' 'lllllU' 1 1 USB I1.U20 oTII 

c) 2S" state wave functions 




'wiiir 1 1 tars* 



.020 0.025 0.030 i 



(d) 25*o wave function differences 



Figure 2: Plots on the left are NQA momentum space wave functions for certain electronic and 
muonic hydrogen states. Plots on the right are differences between NQA wave functions and 
Dirac-Coulomb wave functions for certain electronic and muonic hydrogen states. Thick red lines 
represent electronic hydrogen wave functions or differences and thin blue lines represent muonic 
hydrogen wave functions or differences. 



The full coupled NQA equations are symmetric under mi f-> 771-2 ■ We used an ap- 
proximation to get the final form of the equations used in these numerical calculations. 
This approximation obscures the mass interchange symmetry, but it should still be present 
to some degree. To check this, we interchanged masses and calculated a few of the eigen- 
values for electronic hydrogen, where the mass interchange creates more of a drastic change 
to the equations than in muonic hydrogen. We recovered the same eigenvalues shown in 
the tables up to 1 or 2 sigmas. 

We also found evidence that our precision is not high enough for the final terms in Eq. 



(59) and Eq. (60) to have a significant effect on the eigenvalues. We calculated ground 
state and n = 1 eigenvalues without these terms and the values were not appreciably 
different. This is another motivation for improving our precision in the future. 
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In addition to eigenvalues, we compared our wave functions with Dirac equation 
solutions. Specifically, we compared the momentum space radial wave function of the 
upper components of the Dirac equation solutions with our function Ql{p)- Once again we 
make comparisons for both muonic and electronic hydrogen. Comparisons for two states 
are shown in figure [2] Plots on the left hand side are the NQA wave functions as a function 
of p/mi for three electronic and muonic hydrogen states. Plots on the right hand side are 
differences between the NQA wave functions and Dirac-Coulomb wave functions for the 
same electronic and muonic hydrogen states. We plotted the muonic and electronic wave 
functions and differences for each state on the same set of axes in order to make direct 
comparisons between electronic and muonic hydrogen. As is the case for Dirac-Coulomb 
wave functions, NQA muonic wave functions are larger for smaller momenta and smaller 
for large momenta than their electronic counterparts. It is clear from the right plots that 
muonic hydrogen wave functions are more consistent with Dirac-Coulomb wave functions. 
The NQA solutions are less than the Dirac wave functions for small momenta and greater 
for larger momenta. 

The differences shown on the right in figure [2] are small compared to the size of the 
wave functions themselves, therefore we can conclude that we have found wave functions 
that are fairly consistent with Dirac wave functions, as well as eigenvalues that are all 
within 2 sigmas except for the 1<Sq state. It is worth noting that we introduced the heavy 
particle mass independently in the NQA equations, yet the results compare nicely with 
Dirac equation results with the reduced mass. 



9 Framework for higher order contributions 

In this section, we show how to calculate higher order contributions to the bound state 
energy, specifically Lamb shift terms. The actual calculation of these terms is beyond 
the scope of this paper, but we will explain how to extend our formalism to include the 
corrections. The usual diagrams used to calculate the Lamb shift perturbatively are shown 
in figure 3. The energy contributions of these diagrams are typically calculated with 



respect to zeroth order wave functions of solutions to equations such as Eq. (31). We 
intend to calculate the analogs of these diagrams within the framework of the NQA. We 
can incorporate these terms in our integral equations and use the numerical techniques 
described above to find a less perturbative solution. 

Three examples of NQA Lamb shift diagrams are shown in figure 4. The external off- 
shell line at the lower left of each diagram is assumed to be the electron. These terms will be 
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Anomalous 
magnetic moment 



Vacuum polarization 



Mass re normalization 



Figure 3: Lamb shift diagrams 






(a) 



(b) 



(c) 



Figure 4: Examples of NQA Lamb shift diagrams 



added to the right hand side of Eq. (24). In any diagram, every loop that does not contain 



an NQA amplitude must have one on-shell line. With this rule, there are 2 permutations of 
the vacuum polarization diagram, 3 of the anomalous magnetic moment diagram, and 2 of 
the mass renormalization diagram. The diagrams shown in figure 5 involve the amplitude 
/ e , but there are similar diagrams involving f p where the line directly to the lower left of 
the circle is the on-shell line and the line directly to the circle's lower right is off-shell. This 
means that there are a total of 4 diagrams for vacuum polarization, 6 for the anomalous 
magnetic moment, and 4 for mass renormalization. We also must find similar diagrams to 



add to the right hand side of Eq. (25 ) where the external off-shell line is the proton. 



The NQA seems to be more complicated than the standard procedure where there is 
only 1 Feynman diagram for each Lamb shift contribution, but diagrams of the same type 
are very similar and it is not necessary to calculate each diagram explicitly. Additionally, 
the mass shell delta functions that appear in the NQA simplify calculations. 
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10 Summary and future work 

We used the N-Quantum Approach (NQA) in one-loop order to calculate the energy levels 
and bound-state amplitudes of ordinary and muonic hydrogen and of positronium. We 
can treat other two-body systems, such as the (e/x) system and the (fj,p) system, in an 
analogous way. We used the NQA systematically to find a relation between wave functions 
with the light particle off-shell and the heavy particle off-shell and to find normalization 
conditions for our amplitudes. 

We will use perturbation theory to add corrections to the electron-photon vertex and 
to the photon propagator to include the terms that lead to the Lamb shift. 

In future work we will derive integral equations that include higher-order corrections. 
We will solve these equations numerically without using perturbation theory in order to 
include some of the energy correction terms that are usually calculated perturbatively. We 
will compute the Lamb shift and the hyperfine structure of both electronic and muonic 
hydrogen. We plan to carry these calculations to sufficient order to compare our results 
with the usual methods to see if our methods resolve the muonic hydrogen anomaly. 
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